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ABSTRACT 

The 7r -theorem of dimensional analysis has been used to obtain similarity para- 
meters characterizing the interaction of a charged -particle gas with a magnetic field. 
These parameters are used to derive scaling relations for simulating ion -containment 
systems with electron-containment systems (electron analogs). The scaling conditions 
are used to correlate the results of an electron-containment experiment to a similar ion- 
containment experiment, both reported in the literature. The analysis is extended to 
include a discussion of the possibility of collisionless plasma simulation. Plasma - 
containment simulation is considered from the viewpoint of the stability of static sys- 
tems. The theoretical possibility of flute instability simulation is discussed. 
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ANALYSIS OF ELECTRON ANALOG SIMULATION OF ION 
CONTAINMENT IN MAGNETIC FIELDS 
by Nelson L. Milder and Harold R. Kaufman 
Lewis Research Center 

SUMMARY 

A physical model of a one -component charged fluid was used to develop similarity 
parameters. These parameters, obtained by means of the ir -theorem of dimensional 
analysis, are assumed to characterize the macroscopic motion and containment of the 
fluid in a magnetic field. The similarity described is between an ion -containment system 
and an electron analog of such a system. 

Scaling relations are derived from the similarity parameters. They are used to 
scale ion-containment systems to their electron analogs. These relations are expressed 
in terms of the electron to ion mass ratio m e /m. and the electron to ion experiment 
geometric length ratio i Ji . . 

Mass scaling has not been explicitly introduced in previous analyses. However, the 
results reported herein are in general agreement with two cases reported in the litera- 
ture. These are adiabatic particle motion in a magnetic field and the motion of a col- 
lisionless (Vlasov) charged-particle gas in a magnetic field. The scaling relations were 
used to correlate the results of an electron-model ion-prototype experiment reported in 
the literature. The scaling relations were found to yield good correlation. This indi- 
cates that, under controlled design simulation, dimensional analysis may well be appli- 
cable to interpretation of experimental results. 

An extension of the one -component results to the description of a two -component 
plasma is discussed. The similarity described is between an ion plasma (ion kinetic 
energy much greater than electron kinetic energy) and its electron analog (electron kine- 
tic energy much greater than ion kinetic energy). Here, the simulation of a hot colli- 
sionless ion plasma by means of a hot collisionless electron plasma is examined. It is 
demonstrated that simulation is possible for static systems. The stability of static col- 
lisionless plasmas is considered, and it is shown that simulation of at least one type of 
instability associated with such systems - the flute instability - is theoreticaUy possible. 
This result is in agreement with recent experiments reported in the literature. 



INTRODUCTION 


This paper shows how dimensional analysis can be used in interpreting electron- 
containment experiments as analogs of ion -containment experiments. Such analyses can 
aid our understanding of charged-particle containment phenomena. It extends the appli- 
cation of limited experimental results to physically different systems. 

Electrons have been used for some time to simulate phenomena associated with ion 
motion in magnetic fields. The lower electron inertia and lower magnetic -field strength 
requirements lead to simplifications in the design of the experiments and to substantial 
reduction in costs. As early as 1908, Birkeland studied charged -par tide interaction 
with the Earth’s magnetic field by injecting electrons toward a magnetized sphere (ref. 1). 
More recently, the effectiveness of various magnetic -trap configurations has been studied 
by means of electron injection into model field configurations (refs. 2 to 4). Interest in 
containing high -temperature ion plasmas for production of thermonuclear fusion power 
has led to experiments such as the Table Top IV (ref. 5). In the experiment described in 
reference 5, a hot electron plasma (electron temperature much larger than ion tempera- 
ture) is trapped in an open-ended magnetic well. 

Similarity parameters and scaling relations for plasma devices have been obtained 
by Janes (ref. 6) using the kinetic and Maxwell equations. Schindler (ref. 7) has used 
a similar approach to obtain similarity laws for the interaction of the solar wind with the 
magnetosphere. In the present report, the techniques of the n -theorem of dimensional 
analysis (refs. 8 and 9) are used to establish parameters of similarity between experi- 
ments in which electrons are used to simulate ion motion and the actual ion experiments. 
Such electron experiments can be referred to as electron analogs. 

In the 77 -theorem approach, theoretical and experimental knowledge is used to 
select the pertinent variables. The 77 -theorem, together with these pertinent variables, 
permits derivation of dimensionless variables. These dimensionless variables can then 
be used to simulate an ion device with a much more tractable electron device. This 
technique is advantageous because (1) understanding an ion-plasma device can be attained 
without either the construction of the ion device or an explicit mathematical solution of 
its operation and (2) the performance of related systems (ion and electron) is readily 
correlated, aiding the systematic understanding of such related devices. 

Herein, simple physical models are used to derive similarity parameters and 
scaling relations. The useful properties of these models are that (1) the dimensionless 
parameters can be checked with known equations (In this way the analog concept can be 
verified.), (2) the models can be used to interpret and correlate the results of charged- 
particle containment experiments reported in the literature (e. g. , ref. 10), and (3) the 
models can be used to demonstrate the possibility of constructing plasma analogs in a 
limited sense. 
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A physical model of a one -component charged-particle fluid is presented and analysed. 
The model is used to obtain a set of physical variables from which appropriate dimension- 
less parameters are derived by means of the 77 -theorem of dimensional analysis. These 
parameters are used to establish a set of similarity conditions applicable to designing 
electron analogs of ion trapping experiments. The dimensionless, or similarity, para- 
meters are interpreted in terms of quantities characterizing physical systems. The re- 
sultant similarity conditions are applied to three one -component physical systems. The 
analysis is then used to correlate the data of similar electron- and ion-containment ex- 
periments reported in the literature. Finally, the possibility of using dimensional 
analysis techniques in plasma simulation is discussed. Although explicit application of 
the technique to the problem of plasma containment is not included in this report, the 
possible simulation of the flute instability is considered. 


PHYSICAL MODEL 
Single-Species Charged Gas 

Consider the nonrelativistic steady-state motion of a gas of charged particles in an 
arbitrary magnetic field. The set of equations governing the motion are the set of Max- 
well equations, 


V X E = 0 

(1) 

V. £ = ^ 

(2) 

e o 


O 

II 

1PQ 

> 

(3) 

V x B = M 0 nqv 

(4) 


and the momentum equation (neglecting gravitational body forces). 

nm(v • V)v = nq(E + vxB) - Vp (5) 


(All symbols are given in appendix A; SI units (rationalized MKS system) are used 

throughout this report. ) The kinetic stress tensor reduces to a scalar pressure 
2 

p = 1/3 nmu when collisions are sufficient to produce a Maxwellian particle distribution 


3 



II I Mill Hill IBM II II I i II II mam i i 


with mean random speed much greater than the drift speed (|u| » | v|). The total parti- 
cle velocity V is given by the sum of the random velocity and the net drift velocity 


V = u + v 


( 6 ) 


Here, |u| = (u*0 


( 9 ^ 1/2 

is used for convenience, since the small constant factor discrepancy 
between rms and mean speed is unimportant in present dimensional considerations. It is 
now possible to use equations (1) to (5) to select the pertinent variables that describe the 
physical system. These variables and their dimensions are as follows (The fundamental 
dimensions are mass M, length L, time T, and charge Q.): 


u random or rms particle speed, L/T 
v drift speed of gas, L/T 

B magnetic-flux density, M/QT 

l characteristic length, L 

m mass of charged particle, M 

3 

n density of charged particles, 1/L 

q electric charge, Q 

2 2 3 

€q free space permittivity, T Q /ML 

o 

q free space permeability, ML/Q 

Here the quantity i is used to represent all geometric lengths. The justification for this 
will be considered when conditions of similarity are discussed. The electric field is 
assumed to arise solely from the charge distribution and is therefore not included as a 
separate parameter. 


Application of 7rTheorem 

The preceding set of nine quantities describe the physical system by means of 
functional relations of the form 


<p(u,v,B,z,n,m,q, e Q , ju 0 ) = 0 (7) 

There are four fundamental dimensions in these nine quantities so that, according to 
the jr -theorem, five dimensionless parameters can be formed. Functionally, one has 
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f(7r r n 2 , it 3 , i r 4 , ffg) = 0 


( 8 ) 


The techniques of the n -theorem (refs. 8 and 9) are used to obtain the five dimensionless 
parameters. Making the choice of m, q, and /Uq for the primary quantities, the 
equations to be solved are 


a l D 1 C 1 1 e l 
TTj = u m q e 0 n Q 

(9) 

a 2 b 2 C 2 d 2 e 2 
^2~ B m q € 0 ^0 

(10) 

a„ b„ c„ d„ e„ 
= l m d q 'Vq 

(ID 

a 4 b 4 c 4 ^4 e 4 

?r 4 = n m q e Q 

(12) 

a 5 b 5 n C 5 d 5 e 5 
tt 5 = v m q e Q 

(13) 


The method of solution is to substitute the dimensions of the variables shown and to de- 
termine the values for all the exponents such that the right sides of relations (9) to (13) 
are dimensionless (refs. 8 and 9). In this manner, the resulting dimensionless para- 
meters are 


tt 1 = e 0 p 0 u 


"2 = 




m 


_ _ mZ 
77 3 ^ 


^0 q 


ff 4 = 


3 6 
M 0 q n 


m 


(14) 

(15) 

(16) 
(17) 
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(18) 


*5 = e 0^0 v2 

These parameters may be used in the form shown, or they may be rearranged to obtain 
new dimensionless parameters. The only limitations on the use of new parameters for 
the problem under consideration is that they be five in number and that they should be 
expressible as separate and distinct combinations of ^ 3 * ^ 4 , and The 

separate and distinct combinations chosen for this report, in order to facilitate physical 
interpretation, are 


1/2 


P l = 


P„ = 


"5/ 

rl/2 


u 

V 


mu 


2 w l/2„ Biq 


n 2 n 3 


(19) 


( 20 ) 


P 3 = 


* 1*4 M 0 mnu 


( 21 ) 


ft a 


P 4 = 


P 5 = 


2 

2 2 4 

*1 

e 0 m u 

2 2 

2 2 2 

*3*4 

Mq 1 n <1 


r 2 


e 0 B ' 


( 22 ) 


(23) 


Physical Interpretation 


An understanding of the physical significance of the dimensionless parameters 
(eqs. (19) to (23)) is helpful in evaluating their uses and limitations. The first parameter 
Pj is the ratio between the drift and rms speeds. 


P 


1 " 


(Drift speed) 
(rms speed) 
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The grouping mu/Bq in the second dimensionless parameter (eq. (20)) is effec- 
tively (to within a constant of order one) the cyclotron radius of the charged particle. 
The remaining variable in Pg is the characteristic length of the device. The dimen- 
sionless parameter Pg can therefore be expressed as the physical ratio 



(Cyclotron radius) 
(Characteristic length) 


which is the well known adiabatic ratio. 

o 

The kinetic pressure of the charged -particle gas can be written as mnu /3. The 
transverse stress of a magnetic field is B /2 ^q, so that the third dimensionless para- 
meter Pg can be expressed as the physical ratio 

„ _ 3 (Kinetic pressure) 

2 (Magnetic stress) 

This parameter (without the numerical factor) is generally denoted as (3 in the litera- 
ture. 

The physical meaning of the fourth dimensionless parameter requires a more ex- 
tended analysis. This analysis is included in appendix B and shows that is a meas- 
ure of the transverse diffusional effect of Coulomb collisions. To within a logarithmic 
factor, P^ is the ratio of characteristic length to Coulomb mean free path. 


P 


r>-/ 

4 ~ 


(Ch aracteristic length) 
(Coulomb mean free path) 


o 

In the fifth dimensionless parameter (eq. (23)) the grouping B /2ju n is the magnetic- 
field stress. The remaining grouping n q l /gq can be interpreted to within a constant 
K (which depends on the spatial shape of the charge distribution) as the electric field 
stress. The ratio of the two groupings of variables thus gives 


P 


5 



(Electric stress) 
(Magnetic stress) 


Similarity Conditions 

This report defines conditions of similarity between an electron analog of an ion- 
containment system. To summarize the physical interpretation of the dimensionless 
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parameters, similarity would be established between an ion device and its electron ana- 
log if five physical ratios were equal for the two devices. These ratios are 


P 


_ (Drift speed) 
1 (rms speed) 



(Cyclotron radius) 
(Characteristic length) 


P 3 


3 (Kinetic pressure) 
2 (Magnetic stress) 



(Characteristic length) 
(Coulomb mean free path) 


P 


5 


K (Electric stress) 
2 (Magnetic stress) 


The ratios of geometric lengths can be included as additional dimensionless para- 
meters in an analysis of this type. The approach used herein for establishing similarity, 
however, is that all geometric lengths must scale by the same factor. For example, if 
the length of an electron analog device is half that of the ion prototype, then all other 
electron analog dimensions should also be one -half of the corresponding ion-prototype 
dimensions. For this reason, only one length scaling ratio need be considered between 
ion and analog devices. 


Application of Similarity Conditions to Single-Species Gas 

The similarity discussed in the preceding section can be used to set conditions for 
designing electron analog experiments. Three possible applications are considered in 
which the similarity conditions can be specialized to yield appropriate scaling relations. 
The procedure for deriving scaling relations is presented in appendix C. In each case 
the applicability of relations (19) to (23) is discussed. 

High -density charged -particle gas trapping in ma gnetic field. - In this case the com- 
plete set of equations (19) to (23) are needed. The appropriate scaling relations are 
(appendix C) 
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(24) 


Z 6 M 3 

n. 


B 

B; 


-= Z 3 M 2 


v u 
6 = 6 _ ^ 

V- U. 

1 1 


— = Z M 


(25) 

(26) 
(27) 


Here M is the ratio of electron mass to ion mass m e /m. « 1/2000. It is apparent that, 
for this case and for reasonable ion model dimensions, analog dimensions become so 
large that experiments are impractical. Physically, the difficulty arises from consider- 
ing only a single -species self -interacting gas of charged particles. For this case, large 
electric stresses can develop as a result of space charge buildup, and the parameter 
Pg of relation (23) is particularly significant. 

Low -densi ty, low -temperature case . - In terms of the physical model discussed 
previously, this case corresponds to the situation in which n — 0 and u — 0. The 
Maxwell equations then are not involved in the calculations and the equation of motion of 
the charged particle becomes 


m — = q(v X B) (28) 

dt 

Studies involving the injection of electrons into magnetic fields have been used to 
evaluate the trapping characteristics of various field configurations (refs. 2 to 4). 

Studies have also been used as a basis for checking the validity of the adiabatic theory 
of particle motion in magnetic fields (refs. 11 to 13). In terms of the present simulation 
analysis, these electron-injection experiments can be related, through scaling, to cor- 
responding ion injection and trapping experiments. This is done by specializing para- 
meters (19) to (23) to this situation. 


P 


1 


0 


P 


3 


= 



(29) 

(30) 
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The remaining parameter Pg corresponds to the scale factor derived by Northrop 
(ref. 14). By setting P 2i “ P 2e and using P li - P le’ the following relation is obtained: 


v e _ m iVe u e 
v i m e B i*i Z u i 


(31) 


where 


q. = Ze 

Thus, for example, one can independently scale mass, length, and magnetic -field strength. 
The velocity scaling will then depend on these three scaling factors. 

Low-density, collisionless case . - This case corresponds to the situation in which 
the Coulomb mean free path exceeds all characteristic lengths of the system. The 
equations governing this case are the set of Maxwell equations and the Vlasov equation. 

The scaling laws have been derived by Schindler (ref. 7) using the complete set of equa- 
tions. Scaling relations can be obtained from the present analysis, again, by specializing 
the similarity parameters (19) to (23). For this case = 0 and only Pj, Pg, Pg, 
and Pg need be kept constant for scaling. The scaling relations are 



^ = z 2 ml- 2 


ZML" 1 


(32) 

(33) 

(34) 


These results are in agreement with those of reference 7. 


CORRELATION WITH EXPERIMENT 

It is of interest to attempt to correlate this similarity analysis with experiment. The 
experiment must permit the correlation of data obtained from a scaled electron analog of 
an ion -injection -trapping machine to the data obtained from the actual ion prototype. 
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Such an experiment has been reported in reference 10. Here, an electron-injection ex- 
periment was studied preliminary to the performance of an ion -injection experiment. 

Both experiments consisted of the injection of charged particles parallel to the axis of a 
cusped magnetic -field configuration and measuring the trapping characteristics of the 
configuration. 

The electron model was geometrically scaled by a factor of two over the ion proto- 
type. A 185-eV electron beam was injected into cusped fields with maxima at about 

-4 -4 -4 

90x10 and 135x10 tesla. For the case of 90x10 tesla, the nonadiabatic confinement 

time (time required for an electron to escape as a result of magnetic-field gradients only) 
was found to be 0. 245 microsecond. Prolonged containment was attained by electrostatic- 
ally closing the cusp leaks, so that electron-containment times for B m „_ = 135x10"^ tes- 

nicix 

la ranging from about 1.43 to 5 microseconds were realized, depending on background 
gas pressure. The electrostatic field was designed so as not to penetrate into the trap- 
ping region. 

The ion prototype experiment used proton injection into similar cusped fields with 
maxima ranging from 1. 3 to 1. 7 tesla. Beam injection energy was 17 keV. The energy 
distribution of trapped ions was not determined; however, estimates of the mean energy 
were of the order of several keV. Containment times depended on background pressure 
and ranged from 12 to about 60 microseconds. Extrapolation of the data to zero back- 
ground pressure yielded a nonadiabatic trapping time of about 4 microseconds. 

To examine these results within the framework of the present analysis, it is neces- 
sary to use the similarity parameters (19) to (23). Intuitively, one might initially sus- 
pect that Coulomb collisions in this experimental environment might be negligible (i. e. , 

P^ = 0) since charged-particle densities are relatively low. However, the scaling rela- 
tions derivable from the condition P^ » 0 do not result in correlation between the elec- 
tron and ion experiments. This lack of correlation is consistent with the fact that Cou- 
lomb collisions, although a small effect, can be the dominant source of charged -particle 
diffusion. In order to correlate these experiments by means of the present analysis, it 
is necessary to assume that the parameter Pg vanishes. This parameter is the ratio 
between electric- and magnetic -field stresses, and can be shown to be essentially the 

p 

ratio (v/c) . Thus requiring Pg to vanish, in effect, implies that the condition v « c 
holds, which is the case for most plasma experiments. 

The scaling relations for this case are obtained by a procedure analogous to that 
given in appendix C. 


5- z-W 


(35) 
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(36) 


V Z-J/2 m -3/4 l 5/4 


A. ZM-V2 l 1/2 


(37) 


In steady state, time scales as £ /v, so that 


1 = Z-1/2 m -1/4 l -5/4 


(38) 


r-1 


For values of L = 1/2 and M = 1836, the scaling relations become (Z = 1) 


n i 

J-= 459 
n 


— = 30 


% 

B, 


= 118 


-1= 15.6 
*e 

where <£ and ^ are ion and electron kinetic energies, respectively. Table I compares 
calculations using these relations with the data of reference 10. The data obtained from 
the electron analog experiment were used to calculate corresponding ion prototype values. 
The calculated values are then compared with measurements reported from the ion injec- 
tion experiment. Particle containment times were measured in reference 10 for both 
the electron- and ion-trapping experiments. The nonadiabatic escape time given for the 
ion -trapping experiment was obtained by extrapolating the experimental containment time 
curve (ref. 10) to zero background pressure. 


12 



TABLE I. - COMPARISON OF CALCULATION TO EXPERIMENTS OF REFERENCE 10 


Item 

Containment 

number 

density, 

m" 3 

Mean particle 
energy, 
eV 

j 

Magnetic flux 
density, 

T 

Containment times 

Electron analog 
(experimental 
data, ref, 10) 

' >10 12 

185 

90X10' 4 to 135X10' 4 

Measured with leaks 
open (nonadiabatic 
escape), 

0. 245 \isec, 

Measured with leaks 
closed (diffusional 
escape), 

1. 43 to 5 psec 

Ion prototype 
(experimental 
data, ref. 10) 

5xl0 14 to 5xl0 15 

>10 3 to 17X10 3 

1.3 to 1.7 

Extrapolated to 
zero pressure 
(nonadiabatic escape), 
4 jusec 

Measured as function 
of pressure (differen- 
tial escape), 

20 to 50 (xsec 

Ion prototype (cal- 
culated from elec- 
tron analog data) , 

>4. 6xl0 14 

5.6X10 3 

1.1 to 1.6 

Nonadiabatic escape, 
3.8 /i sec 

Diffusional escape, 
22 to 78 iisec 


Experimentally, the containment time increased with background pressure. Within 
the framework of the present analysis, this time varies as I /v. Momentum -transfer 
collisions between charged particles and background particles would effectively decrease 
v, thus increasing t. Such collisions result in the diffusional losses observed in the 
experiments. 

The good agreement between calculated and measured ion -prototype variables would 
seem to indicate that correlation between an electron analog of an ion-trapping experi- 
ment using dimensional analysis techniques is indeed possible. 


PLASMA SIMULATION 

So far we have considered systems containing positively or negatively charged par- 
ticles only. It may be even more interesting to determine whether these similarity 
parameters can be used to scale plasma systems; if so, what are the limitations of such 
a simulation? 

For purposes of discussion, one can distinguish between an ion plasma and an elec- 
tron plasma. The ion plasma is characterized by ion energy much greater than electron 
energy; that is, electrons serve as space -charge neutralizing background particles. 

For an electron plasma, the situation is reversed (electron energy is much greater than 
ion energy). Herein, the conditions for which an electron plasma can be regarded as the 
analog of an ion plasma will be determined. 



In the following analysis, two -component plasmas are considered. The two species 
are regarded as interacting fluids, and the appropriate hydrodynamic equations are in- 
vestigated. A complete set of equations describing these systems are as follows: 


Momentum equations: 

/Sv 


n m „ | 
as as 


«£ + v 


at. 


as 


Vas - n as<W g + Vas * 5 s) “ V«s “ 5 a/3, s < 39 ) 


Continuity equations: 


Maxwell equations: 




V S " E S ^ n Q!S^Q!S 


: 0 a 


V ■ B = 0 
s s 


v s xE s = 


0t„ 


V.XB = 


M 0 


V - 1 3E s 

a n ^s v as + ^ ^ 


(40) 

(41) 

(42) 

(43) 

(44) 


where 



D s = e 0 E s 


Equations of state: 


n = i n m u 3 _ 
H as 3 as as as 


(45) 


In these equations the subscripts a and j3 denote either ions or electrons (a / $) and 
s denotes the system (ion plasma, s = 1; electron plasma, s = 2). 
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Now, consider the following problem. Under what conditions can the previously de- 
rived similarity parameters be used to transform system s = 1 equations into system 
s = 2 equations, leaving the form of the resulting equations invariant? The analysis 
will again be restricted to steady-state systems. This eliminates all 9/3t terms in the 
equations. Two sets of physical variables must be considered for each plasma system. 

If a denotes the primary species and j3 the background species, then the two-species 
variables are n _, n. . v . v\ . u u„ , m . m. . (P 0 is considered sepa- 
rately. ) In order to bring the problem into the realm of the one -component similarity pa- 
rameters (19) to (23), equations (39) to (43) must be appreciably simplified. To begin with, 

consider the term P „ . This term represents the momentum transferred by colli- 

Gfp, ® 

sions per unit time between the primary charged species a and the background species 
/3 in system s. From conservation of momentum, we have that P^ g = - P^ g . In 
the present case it will be assumed that collisional mean free paths are sufficiently long 
to permit neglect of these momentum -transfer terms. Thus the analysis which follows 
holds for collisionless plasmas. Additional assumptions required are as follows: 

(1) Electron and ion number densities are equal for both s = 1 and s = 2. 

® W'Vs « m e / m i 

< 3 >%S»V 

(4) Second order terms in the equation of motion for the background species are 
negligible. 

(5) q as = " q /3s ( e * » el ectrons and protons) 

Equations (39) to (45) then reduce to the following set: 


Momentum equations: 


n a,s m a/<*s • Vos = n <M q <*s (i s + 7 as x § s> ‘ Vos 

(46) 

0 = VW 2 s - 1 - % x 5 S ) - Vfte 

(47) 

Continuity equations: 


V s ’ ^ n o!S V ois^ = V ^ n /3s v /3s^ = ® 

(48) 

Maxwell equations: 


V • E = 0 
s s 

(49) 

v s • 5 s = 0 

(50) 

v s x E s “ 0 

(51) 

V s X 5 S - '“o’Vs'W 1.8 

(52) 
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Equations of state: 


P = — m n u 2 

r A| g _ me /uC / 


as as as 


and 


Ops = l m ss V u Is 


( 53 ) 


From assumption (2), and equation (53) the partial pressures p aS and p^ g are such 
that p ag » p^g. Adding the momentum equations (46) and (47) yields the total momen- 
tum conservation equation 


n m (v • V v ) - (v „ - v _) XB„ - V p o 
s ofS v as s as' ' a s /3s' s s s 


(54) 


The only Maxwell equation involved is 


V s xB s = ^Ws^s - >) 


(55) 


Note that the change of variables p g = Pq , s + P/3g - p aS , n g =n ag = n^, and q g = q^g = 
- q^g has been made. From assumption (3), we have that (v ftS - v^ g ) « v^ g in equa- 
tions (54) and (55). In this manner the complete set of plasma equations can be reduced 
to one involving single -particle variables only, and the similarity parameters (19) to (23) 
are thus applicable to such systems. For nonrelativistic particle speeds, Pg can be 
neglected (see the argument given in the Experimental Correlation section). The scaling 
relations appropriate to this plasma system are as follows: 


m 


It 2 


= Z M L 
n 0 


(56) 


!!i= 

B„ 


(57) 


l v il - v ell 
f e2 - f i2| 


v. 


il _ „1/2 m 1/4 t 1/4 _ u il 


■= Z J 


J M 


e2 


u 


e2 


( 58 ) 


where s = 1 denotes the ion plasma and s = 2, the electron analog, 
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and 



Using the steady -state time scaling t ~ l/v results in 


*1= z 

to 


(59) 


By solving for all the quantities with subscript 1 in terms of subscript 2 quantities 
and substituting into equations (54) and (55) for s = 1 (noting that = LV^), the ap- 

propriate s = 2 equations are obtained. They are identical in form to the s = 1 equa- 
tions. Thus the total momentum, continuity, and Maxwell equations for the described 
plasma systems are invariant under the proposed transformations (and assumptions). 


Possible Simulation of Plasma-Containment Systems 

Of great interest would be a rigorous demonstration of the possibility of using hot 
electron plasma containment to simulate hot ion-containment systems. Difficulties are 
encountered from the fact that many observed and predicted instabilities associated with 
such systems result from locally, or microscopically, developed growth mechanisms 
which cannot be simply correlated with macroscopic variables. For the dimensional 
analysis used in this report, the complexity of the relation between instability growth 
mechanisms and physical variables need not be known exactly. However, it is essential 
to know all the variables on which these mechanisms depend in order to obtain appro- 
priate similarity parameters. There does exist a class of instabilities for which the 
nine quantities of equation (7) are applicable, and thus the scaling parameters derived 
herein can be used to describe their simulation. 

Simulation of Flute Instability 

One type of hydromagnetic instability which is known to depend, in part, on the 
magnetic -field configuration is the interchange, or flute, instability. This instability 
has been extensively considered in the literature (refs. 15 to 17). An energy principle 


17 



has been used (ref. 17) to derive the following stability criterion for open-ended mirror 
systems. 



(60) 


Here, the integration is performed along a magnetic -field line. The radius of curvature 
R of a field line is defined as positive at the mirrors. The radius r is measured from 
the symmetry axis. Consider this integral in an analog experiment. The scaling rela- 
tion (57) can be used to derive the following relation between model and prototype 


f df e (L f 2 f df i 

J r R Bl ^ J r.R.E 

e e e 


(61) 


r i R i B i 


Since the scale factor (L/M)^^ is always positive, it is evident that, if the left side of 
equation (61) is positive, the integral on the right side will also be positive. Thus, si- 
mulation of flute stabilization should be theoretically possible for mirror systems. This 
conclusion is in qualitative agreement with results obtained from recent Table Top IV 
experiments (ref. 5). Therein a linear quadrupole magnetic field superimposed on an 
axial field stabilized a hot electron plasma against the flute. Although scaling was not 
rigorously satisfied, the authors concluded that their results (including a measured 
instability) may have been the electron analog of hot ion plasma, quadrupole -stabilized 
systems such as DCX, OBRA, PHOENIX, and ALICE. 


CONCLUSIONS 

A physical model of a one -component charged fluid was used to develop similarity 
parameters. These parameters, obtained by means of the it -theorem of dimensional 
analysis, are assumed to characterize the macroscopic motion and containment of the 
fluid in a magnetic field. The similarity described is between an ion-containment sys- 
tem and an electron analog of such a system. 

Scaling relations are derived from the similarity parameters. They are used to 
scale ion -containment systems to their electron analogs. These relations are expressed 
in terms of the electron to ion mass ratio m e /m. and the electron to ion experiment 
geometric length ratio i /i 
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Mass scaling (e. g. , replacing protons by electrons) has not been explicitly intro- 
duced in previous analyses. However, the results reported herein are in general agree- 
ment with two cases reported in the literature. These are adiabatic particle motion in 
a magnetic field and the motion of a collisionless (Vlasov) charged-particle gas in a 
magnetic field. The scaling relations were used to correlate the results of an electron 
model -ion prototype experiment reported in the literature. The scaling relations were 
found to yield good correlation. This indicates that, under controlled design simulation, 
dimensional analysis may well be applicable to interpretation of experimental results. 

An extension of the one-component results to the description of a two-component 
plasma is discussed. The similarity described is between an ion plasma (ion kinetic 
energy much greater than electron kinetic energy) and its electron analog (electron 
kinetic energy much greater than ion kinetic energy). Here, the simulation of a hot col- 
lisionless ion plasma by means of a hot collisionless electron plasma is examined. It 
is demonstrated that simulation is possible for static systems. 

The stability of static collisionless plasmas in a magnetic field is also considered, 
and it is shown that simulation of one type of instability associated with such systems - 
the flute instability - is theoretically possible. This result is in agreement with recent 
experiments reported in the literature. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 27, 1969 
120-26-02-01-22. 
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APPENDIX A 
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SYMBOLS 

exponents in relations (9) to (13) 
magnetic -flux density 
impact parameter in appendix B 
exponents in relations (9) to (13) 
distance of closest approach, q / 43 T€Qmv 
speed of light 

exponents in relations (9) to (13) 
displacement current 
exponents in relations (9) to (13) 
electric field 

particle kinetic energy in eq. (37) 
exponents in relations (9) to (13) 
magnetic -field intensity 
geometrical shape factor in Pg 
length ratio 

length measured along field line 
characteristic length of system 
effective cut-off length for Coulomb interactions 
mass ratio 

charged-particle mass 
charged -particle density 
momentum transfer term 

dimensionless parameters defined by relations (19) to (23) 

pressure 

electric charge 

radius of curvature of a field line 


20 



r 

t 

u 

u 

V 
v 

V 

Av l 

(AVj/v ) 2 

6v l 

w 

Z 

e o 

X c 

Mq 

<P 


radius measured from symmetry axis in eq. (60) 
time 

rms speed 

rms velocity 

total particle velocity 

drift speed or particle speed 

drift velocity 

total velocity change averaged over many encounters 

dimensionless quantity proportional to diffusion coefficient (herein referred 
to as a scattering parameter) 

incremental change in velocity normal to initial velocity direction resulting 
from a single binary encounter 

mean relative speed of incident test particle in appendix B 

electronic charge number 

free -space permittivity 

2 2 2 2 

Coulomb mean free path, 47reQ(mv ) /q In A 
free -space permeability 

dimensionless parameters defined by relations (14) to (18) 
generalized function defined by eq. (7) 


Subscripts: 

e electrons in analog 

i ions in simulated system 

s plasma system 

a electrons or ions 

/3 electrons or ions ( a ^ 0) 
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APPENDIX B 


THE PHYSICAL INTERPRETATION OF P 4 

Consider a charged particle with mean relative speed w traversing a gas of like 
particles of uniform density n. The mean kinetic energy of the gas is assumed to be 
much less than the kinetic energy of the incident particle, so that the mean relative speed 
w is very nearly the incident particle speed v. As a result of Coulomb collisions with 
the background particles, the incident particle will suffer a cumulative large angle de- 
flection after traversing a characteristic length l of the gas. The total deflection is 
measured by the quantity (AVj/v) 2 , where Av ± is the total deflection normal to the 
incident direction. Herein the quantity (Av ± /v) 2 is defined as a scattering parameter 
which can be related to a diffusion coefficient. 

For a binary encounter between two charged particles, the incremental change in 
velocity normal to the incident direction is given by 


5v, 


27re 0 mbv 


(Bl) 


where b is the impact parameter. Squaring equation (Bl) yields 



47r 2 e 2 m 2 b 2 v 4 


(B2) 


Denote by n the background density resulting from some space charge distribution es- 
tablished in a time of the order of l/v. Then the scattering parameter is given by 

0 

The volume element dr is that swept out by the test particle as it passes through the 
scatterers of uniform density n. Thus, 

dr = 2nib db 

It is known from energy considerations that the lower limit on the integral must approxi- 
mate the distance of closest approach 
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u 2 
4776Qmv 

The upper limit I ^ denotes an effective cut-off length for Coulomb interactions. Gen- 
erally, this length is taken to be the Debye length (refs. 18 and 19). This is rigorously 
true for a plasma in which thermal equilibrium is maintained. For the present problem, 
in which there are charges of only one sign, there is no cut-off distance short of the 
extent of the charged-particle gas I . It is therefore assumed that Zg — Z in equa- 
tion (B3), so that substituting equation (B2) into (B3) and performing the integration yields, 
for the scattering parameter, 



q 4 nz 
2jre 2 m 2 v 4 


In 



(B4) 


This relation can be written in terms of the dimensionless parameters Pp Pg, P^, 
and P 5 given by relations (19) to (23). Substituting for b Q and multiplying the numera- 
tor and denominator of the logarithmic factor by (n z) 4//2 yields the result 



— — in 

2!7€gm 2 V 4 


|47re 0 mv 2 (nZ 3 ) 

qV/V / 2 


(B5) 


so that in terms of the dimensionless parameters 



(B6) 


Because Pp Pg, P 4 , and Pg are independent quantities, P^ can be altered without 
affecting the ratio Pg/P 2 Pg in the logarithm. Thus the scattering parameter has the 
form 



if C = 4jr(Pg/P 2 Pg) is held constant. 


(B7) 
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The scattering parameter (B7) can be shown to be equal to the ratio of the character- 
istic length l to a mean free path. The mean free path is defined as that corresponding 
to a 90° deflection of the incident electron. The total deflection of the particle after 
traversing a characteristic length l of gas is given by (Av^. The number of colli- 
sions required to deflect the particle through a right angle is the number of deflections 

Av. required to change all the particles' motion into transverse motion. This number is 
1 2 

(v/Av.) . The characteristic time for traversing a distance l is l/v. Thus the mean 
0 

free path for 90 scattering is 



9 2 2 4 

27T€Qin v 


q n In 


(B8) 


It is thus apparent that 



(B9) 


and thus the scatter ing parameter and the ratio i/\ have the same functional depend- 
ence on P^. 
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APPENDIX C 


DERIVATION OF THE SCALING RELATIONS 


In the following calculations, we use the similarity parameters (19) to (23) to obtain 
scaling laws. These scaling laws are applied specifically to electron simulation of ion 
motion in a magnetic field. The subscripts on the physical quantities distinguish between 
electron and ion quantities. Thus m g is the electron mass, and m^ is the ion mass. 

Similarity between an electron analog and an ion prototype is rigorously established 
if 

p e j = p ij i = 1, 2, 3, 4, 5 (Cl) 

Using equations (19) to (23), we obtain from (Cl) 
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e 0 B e 

e 0 B ! 2 


(C2) 


(C3) 


(C4) 

(C5) 


(C6) 


Now, q g = e and q. = Ze. Also, £q and juq are the same for both electron and ion 
systems. Equations (C2) to (C6) then yield 
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B? 


From equation (C9), 



From equation (C8), 



Substitute equation (C12) into equation (C13) to obtain 



(C7) 


(C8) 


(C9) 


(CIO) 


(Cll) 


(C12) 


(C13) 


(C14) 
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From equations (CIO) and (C14), 



Substituting equation (C15) into equation (C12) leads to the result 


(C15) 



u i 


Substituting equations (C14) and (C16) into equation (Cll) yields, finally, 


(C16) 


1= Z 2 ^ 
e m i 


(C17) 


Equations (C14) to (C17) constitute the complete set of scaling laws satisfying the simi- 
larity conditions (Cl). If the mass ratio is defined by the scaling factor M, 

then equations (C14) to (C17) can be written as follows: 


Z 6 M 3 

n. 


?£= Z 3 M 2 
B i 


u 

e _ ^ _ g 

u. v 

1 v i 


— = z 2 m 


(CIS) 


(C19) 


(C20) 


(C21) 


For those cases in which the assumption of space charge neutralization is valid, similar- 
ity conditions Pg. = Pg e can be neglected. Proceeding in a manner analogous to that 
used with relations (C7) to (C9), one obtains the following scaling relations 
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(C22) 


(C23) 


(C24) 


Introducing the scaling factors M and L results in relations (C22) to (C24) becoming 
relations (35) to (37) of the text (Equation (37) is obtained by multiplying equation (C23) 
by M and inverting). 
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